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Abstract 

An analytical expression for the first-order density matrix of a charged, two-dimensional, har- 
monically confined quantum gas, in the presence of a constant magnetic field is derived. In contrast 
to previous results available in the literature, our expressions are exact for any temperature and 
magnetic field strength. We also present a novel factorization of the Bloch density matrix in the 
form of a simple product with a clean separation of the zero-field and field-dependent parts. This 
factorization provides an alternative way of analytically investigating the effects of the magnetic 
field on the system, and also permits the extension of our analysis to other dimensions, and/or 
anisotropic confinement. 
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I. INTRODUCTION 



Theoretical investigations of harmonically trapped ideal Fermi gases have seen a renewed 
interest in recent years owing to the remarkable experimental advances made in the area of 
trapped, ultracold atoms. 1 Indeed, theorists now have an experimental realization of what is 
close to being an ideal, inhomogeneous, quantum many-body Fermi system. Sophisticated 
magneto-optical traps now allow for the possibility of "tuning" the dimensionality of these 
gases from three dimensions (3D) to quasi-2D or quasi-lD. Thus, studying the properties of 
essentially ideal, lower-dimensional many-body Fermi systems is now firmly in the realm of 
experimental fact, and not simply a matter of academic interest. Furthermore, analytical 
results for these systems can be of great use in the density-functional theory (DFT) of inho- 
mogeneous Fermi systems, whereby one can bypass the numerically expensive one-particle 
Schrodinger equations. 2 

The ideal charged Bose gas (CBG) is the Bose analog of a charged Fermi system. This 
model consists of a gas of spinless, charged bosons, coupled to an external, homogeneous 
magnetic field, and was first investigated in 3D by Osborne, 3 Kosevitch, 4 and later by 
Schafroth, 5 The uniform 2D CBG has also been analytically studied quite extensively in the 
literature in light of its possible connection to the theory of high-T c in the cuprates . 6 ^ 8 ^ 10 ' 11 
To date, no detailed analytical analysis has been performed for the inhomogeneous case. 
Since the confined 2D CBG is no longer forbidden from undergoing a Bose-Einstein con- 
densation (BEC) transition at low temperatures (i.e., the Bogoliubov 1/k 2 theorem is no 
longer applicable),— an exact analytical investigation of the thermodynamic and magnetic 
properties (e.g., the Meissner-Ochsenfield (M-O) effect) of the inhomogeneous system would 
be of great interest. 

The fundamental quantity from which the thermodynamic and magnetic properties of the 
ideal quantum gases are derived is the first-order density matrix (FDM), p(r, r'). However, it 
is highly non-trivial to obtain an exact expression for the FDM (even at zero temperature) 
for all but the simplest of cases, viz., the homogeneous ideal charged quantum gas. The 
introduction of a magnetic field further complicates the problem, and it is only relatively 
recently that an exact analytical expression for the zero temperature FDM of a uniform 
Fermi system coupled to a homogeneous magnetic field has become available.— Extensions 
of these results (i.e., to include the case of the CBG and finite temperatures), have only 
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been given in the last few years.— For the inhomogeneous ideal charged quantum gas, even 
the field-free case at zero temperature is difficult. Indeed, exact results for non- uniform 
systems at zer o 14 i 15 i 16 i 17 anc l finite temperatur o 18 i 19 i 20 are limited to the case of harmonic 
confinement. To our knowledge, closed form, exact results for p(r, r') for an ideal charged 
Fermi or Bose gas under general confinement, finite temperatures, and arbitrary magnetic 
field strength, are not available. 

The purpose of the present work is to help fill in this gap by providing an analytical ex- 
pression for p(r, r') , generalized to treat exactly the presence of a uniform external magnetic 
field and confining potential. Our focus will be on providing results for the 2D harmonically 
confined quantum gas, although the general approach of our analysis does allow for an exten- 
sion to other dimensions, and anisotropic traps, should the need arise.— The exact results of 
this paper should prove useful in the areas of current-density-functional theory, (CDFT ) 22 ' 23 
which is a rigorous extension of DFT to inhomogeneous systems immersed in an external 
magnetic field, and for the analytical investigation of the magnetic and thermodynamic 
properties of the CBG in the case of nonuniform systems. 

The rest of our paper is organized as follows. In the next section, we introduce the central 
theoretical tool used in our analysis, viz., the Bloch density matrix (BDM). In, Section III, 
we provide a derivation of the inverse Laplace transform of the BDM, which leads directly to 
the exact FDM for a Fermi or Bose gas at any finite temperature and magnetic field strength. 
Section IV summarizes our main results and offers a discussion of how they may be applied 
in the context of CDFT and the inhomogeneous 2D CBG. 

II. THE BLOCH DENSITY MATRIX 

The central theoretical tool used in our analysis is the zero temperature BDM, Cq(t, r'; /3), 
which is related to the FDM through an inverse Laplace transform.— One of the key reasons 
for working with the BDM is that one does not require explicit knowledge of the one-particle 
wave functions of the associated trapping potential. In addition, the zero temperature BDM 
is independent of the quantum statistics of the system, thereby allowing for an extremely 
robust approach for treating either the Fermi or Bose gas. Since a detailed discussion of the 
BDM has already been given in our previous work, we will only present here the essential 
formalism required for a self-contained statement of the problem, and refer the reader to 
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Refs. ll0lll4Ul8lll9l for additional details. 
The zero temperature BDM is defined by 



C (r,r';/3) = £ # (r')^(r) exp(-/fei) , (1) 

alii 

where the ^j's and the one-particle energies e, are solutions of the Schrodinger equation. The 
constant (3 above is to be interpreted as a mathematical variable which in general is taken 
to be complex, and not the inverse temperature 1/ksT. The BDM satisfies the so-called 
Bloch equation 

HMl{l y, J) = J^01, (2) 

subject to the initial condition 

C (r,r';0) = 5(r-r') . (3) 

In this paper, 

(p - eA/c) 2 1 , , 2 2 . , iS 

is the specific Hamiltonian we work with, where the magnetic field B = V x A, is applied 
along the z-axis, and 

A=(- l -By, l -Bx,^ . (5) 

Note that while Co(r, r'; (3) is gauge dependent, any physical observable is necessarily gauge 
invariant. By choosing a general functional form for Co(r, r'; the solution to Eqs. (2-3), 
with the Hamiltonian (4), can be obtained without having to specify the single-particle wave 
functions or energies. Such a solution has already been obtained by March and Tosi,— which 
we now present in a more explicit form: 

_ . , mUJeff 1 _ ' m "cff sinh(ft^ e /3) ( I I) 

C (r, r ;/3) = — e h ^{^ e «P) Kxy v ' 

2nh smh.(hiu e f[(3) 

x e -[(,-,o 2 +fa-y0 2 ]^[coth(^ cff/3 ) + -^t7 ff % 

x e -[(x+x0 2 +(y+y0 2 ]^[c°th(^ cff ffl- B ^^% J ; (()) 



where 



eB k 



Uc= 2^~c' W ° = V m ' u <& = V u o + ^1 . (7) 
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Introducing the center-of-mass and relative coordinates q and s, respectively, 

r + r' , 
q=— , s = r-r', (8) 

allows us to write the BDM as 

C q,s;^ = — — e h sJn ^« 

27TAI smh(ftw e ffp) 



x e 



ft 



x e~ 



For later convenience, we now introduce the following definitions: 

5^ 1 /" \ i id 

A = q 2 + ^-, B=-\—-q 2 \--(q y s x -q x s y ), W = ^' ( 10 ) 

and scale all lengths and energies by £ e fr = \/^/mWrft and hoo e g, respectively. The zero 
temperature BDM can then be written in the more compact form 

C (q,s;/3) = - - exp (-A coth(/3) - B^— - B*-^—\ , (11) 
27rsinh(p) \ sinh(p) smh(fj) J 

where B* denotes complex conjugation. Equation (11) serves as the starting point for the 
rest of our study, but it is worthwhile pointing out that a novel factorization of the BDM 
can be performed. 

First, let us re- write the BDM as follows 

2 s2 \ t_ / o\ ( 2 s2 \ coth(u;/3) 



Co(q,s;/3) = - — 77^ exp t - [q 2 + - coth(/3) + [q 2 - - 



27rsinh(/?) ^ [ V 4 / V 4 / sinh (/?) 

sinh(u;/3) 1 

+ Hgy^ ~ g»gy) . S . (12) 



sinh(/3) 

Making use of the trigonometric identities 

c 05 h(^ -coshM) = _ 2 !^W) 

smh(p) smh(p) 

cosh(/3) + cosh(u;/3) , , sinh 2 (u;/5/2) 

^ 1 ^ = coth (/3 2) + 2 , \ p ' } , 13 

sinh(p) sinh(p) 

in Eq. (12) gives 

Coins; 13) = l^ e h 2 ^m)-icot Hm } eUM) (M) 

27rsinh( / 9) 



5 



where 



rrf ^ n ( 2 s 2 \ smh 2 (tuP/2) .. . sinhM) . ir . 

The "effective potential" U c (q, s; (3) explicitly includes all of the magnetic field dependence, 
and there is a clean separation of the BDM into field free and field dependent parts. In 
particular, setting uj c = in Eq. (15) immediately gives U c = and the BDM, (14), reduces 
to that of a 2D harmonically trapped system . 14 i 18 i 19 This factorization is reminiscent of 



26l which was motivated by the desire to 



the introduction of an effective potential in Ref. 
improve the Thomas-Fermi approximation to potentials which are varying too rapidly in 
some regions of spaced Viewing the magnetic field as an additional ID confining potential 
suggests a similar interpretation in the present context; that is, going beyond uo c = may 
be achieved through the introduction of some effective potential, U c , which encodes the 
magnetic field dependence. Irrespective of this suggestive connection however, Eq. (14) 
here should be viewed as a more direct route to generalizing our results below to other 
dimensions, and allowing for a more transparent analytical investigation of the effects of the 
magnetic field in the weak/high field limits, along with anisotropic confinement, should this 
be desired. 

We are now in a position to see why the BDM provides such a universal scheme to 
investigate either the Fermi or Bose gases. While the zero temperature BDM is independent 
of the quantum statistics, at finite temperature, the BDM for the Fermi system is obtained 
via (k B = 1)24 

C r (q,s;/3) = Co(q,s;/3) . . , (fermions) , (16) 

sm(7rpT) 

whereas for bosons, it is given by 

CV(q, s; 13) = C (q, s; P) ^^^ > (bosons) . (17) 

Therefore, aside from the different temperature dependent factors in Eqs. (16,17), it is clear 
that only the T = BDM is required to study either quantum gas. 
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III. THE FIRST-ORDER DENSITY MATRIX 



A. Fermi gas 

The (spin-averaged) FDM at finite temperature is obtained by a two-sided inverse Laplace 
transform of the finite temperature BDM. The inverse Laplace transform must be two-sided 
to allow for the dual variable to go negative. Specifically, we have^ 

p(q,s;T) = £- 1 



|cV(q,s;/?) 



(18) 



where [i is the chemical potential, which at fixed uj c , is determined by particle number conser- 
vation. As we have discussed before , 14 i 18 i 19 it is very difficult to perform the inverse Laplace 
transform by simply substituting the finite temperature BDM, as given by Eqs. (11,16), into 
Eq. (18). In order to proceed any further analytically, one requires the following identities 

oo 

exp(-Acoth(/3)) = ^L fc (2A)e- A {e~ 2fc/3 - e~ 2(fc+1)/3 } (19) 

k=0 



-P(-S)=g^l J{e--e~ } ,20, 




exp 



n=0 j=0 J ' \ n ~J 

where Li(x) is a Laguerre polynomial. Utilizing these identities in Eq. (11) gives 



/ m 


M 


f n \ 


\ m — i , 




{n-j I 



Z=0 m=0n=0 i=0 j=0 J 
x | e (-2Z-2m-2n+i+j)/3+(j-i)^ _ £ e (-2l-2m-2n+i+j-2)l3+(j-i)up^ 

X |3 e (-2i-2m-2n+j+j-4)/3+(i-i)o;/3 _ e (-2l-2m-2n+i+j-6)p+(j-i)w/3 j ^2) 

Notice that all of the /3 dependence is now contained in the exponential factors and the 
(two-sided) inverse Laplace transform is now tractable. The mathematical details of this 
transform closely follows our earlier work , 1 ^ 1 ^ 1 ? so here we will simply give the final result, 
namely, 

1=0 m=0 n=0 i=0 j=0 k J ' 

x $ k (l,m,n,i,j) , (23) 
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/ m 


II 


C n \ 


\m — i , 




{n-j j 



where all of the temperature dependence is encoded in the Fermi-like function 

1 



(>X p | k+2(l+m+n)-i-j-(J-i)u-iJ, j _|_ ^ 



(fermions) 



(24) 



and the /c-summation is over k = 1,3,5. Equation (23) is the central result of this paper 
and gives the exact FDM for an ideal, harmonically trapped 2D Fermi gas at arbitrary 
temperature and magnetic field strength. Putting s = immediately yields the spatial 
density of the system. It is certainly worthwhile re-emphasizing that all previous results 

4 



2 







— X 

"*>, \\. 












1 


\> * - 
i — __- 







2 4 
r 



FIG. 1: Plot of the spatial density for N = 110 fermions at zero-temperature and various magnetic 
field strengths. The solid curve is for uj = 0, the dashed curve for u = 0.45 and the dotted curve 
for u = 0.70. All lengths and energies have been scaled as discussed in the text 

pertaining to this system can now be obtained from (23) upon taking various limits. For 
example, in the uniform case at zero temperature, Eq. (23) can be shown to reduce to (with 
dimensional constants recovered) 

„( q+ |, q -|) = ^e^m^/W^^e-^^LW^ (^s 2 ) . (25) 

As an illustrative numerical example, we present in Fig. 1, the spatial density for N — 110 
particles at zero-temperature and various magnetic field strengths. It is important to note 
that while the summations at finite temperature in Eq. (23) look somewhat formidable, any 
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practical numerical implementation requires only a relatively small number of terms. Figure 
1, for example, required only a few minutes to plot using a nominally equipped PC running 
a generic flavor of Unix. Consequently, finite temperature effects are readily studied, should 
the need arise.— Note also that the relative ease for which we were able to write down 
p(q, s; T) should not be used to conclude that the calculation is trivial. In particular, it is 
notoriously difficult to treat finite temperature effects exactly in the Fermi gas owing to the 
fact that one cannot express the Fermi distribution function as a convergent power series, 
except at very high temperatures.- 20 Rather, one should view our almost immediate statement 
of the FDM as a testament to the utility of the inverse Laplace transform technique. 



B. Bose gas 

The power of the inverse Laplace transform technique is also apparent if one wishes 
to extend our results to the case of a harmonically confined 2D CBG. Indeed, one can 
immediately write down the final expression for the FDM, with the only changes being a 
change in the sign in front of unity in the denominator of Eq. (24), viz., 

$Jl,rn,n,i, ?') = -, r (bosons) , (26) 

V ' ' ' ,JJ cxp ( k+2(l+m+n)-i-j-{j-i)u-A l V !) V ' 

and the elimination of the factor of two in Eq. (18) (i.e., the bosons are spinless). Thus, 
with no additional work, we also have an exact, closed form expression for the FDM of the 
trapped 2D CBG at arbitrary temperature and magnetic field strength. Of course, setting 
ujq = reproduces the recently obtained finite temperature results corresponding to the 



uniform 2D CBG found in Ref. 
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IV. SUMMARY AND FUTURE WORK 



We have derived an analytical expression for the FDM of an ideal, harmonically trapped 
charged 2D Fermi or Bose gas at finite temperature and arbitrary magnetic field strength. 
To our knowledge, this is the only example where such an exact expression has been obtained 
for an inhomogeneous quantum gas. Aside from their inherent technical merit, we believe 
that our results now open up several other fruitful avenues of investigation, which are outside 
the intended scope of this paper. 
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First it would be interesting to investigate the properties of the 2D Thomas-Fermi ki- 
netic energy functional for the case of finite magnetic field. The motivating factor behind 
this suggestion lies in the remarkable fact that for zero-field, the 2D Thomas-Fermi kinetic 
energy functional leads to the exact quantum mechanical kinetic energy (i.e. without gra- 
dient corrections) when integrated over all space. This non-trivial result was only recently 
discovered by Brack and one of us.— From Eq. (25), it can be shown that the local-density 
approximation (LDA) 2 of the magnetic 2D kinetic energy functional is identical in form to 
the zero-field case, but now the magnetic field is encoded implicitly by the densit y 10 ! 13 Thus, 
determining whether the 2D magnetic-Thomas-Fermi kinetic energy functional is also exact 
(i.e., similar to its zero-field counterpart without gradient corrections) would be very inter- 
esting. It would also be illustrative to study the 2D exchange energy (i.e., suitable for the 
study of parabolically confined quantum dots in a magnetic field) via the exact FDM and 
compare the integrated, and spatial properties against the commonly used LDA of CDFT. 
This type of comparison has already been undertaken for the zero magnetic field case ,— 
and would be an equally worthwhile endeavor for the finite-field case. Furthermore, having 
an exact expression for p(r, r'; T) also allows for the perturbative study of the effects of 
particle-particle interactions, similar to what has already been performed for the uj c = 
case in Ref . |3 

As for the CBG, the most obvious application of our results will be in investigating 
the thermodynamic and magnetic properties of the inhomogeneous system. In contrast to 
the uniform 2D CBG, the trapping potential stabilizes the system to density and phase 
fluctuations and allows for the possibility of a transition to a BEC— We recall here that 
the uniform 2D CBG does exhibit an essentially perfect M-0 effect, in spite of the absence 
of a BEC state.— Thus, analytically studying the connection between the onset of the M-0 
effect and the BEC phase in the trapped system is, in our opinion, an important problem. 
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